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Preface 


The  U.S.  Army  Corps  of  Engineers  Committee  On  Tidal  Hydraulics  held 
its  101st  meeting  in  Philadelphia,  PA,  on  9-11  September  1992,  at  the  invita¬ 
tion  of  LTC  Kenneth  H.  Clow,  Commander  of  the  U.S.  Army  Engineer  Dis¬ 
trict,  Philadelphia.  The  principal  purpose  of  that  meeting  was  to  review  the 
situation  at  Indian  River  Inlet,  Delaware,  and  to  assist  the  District  in  evaluat¬ 
ing  generalized  diannel  scour  there. 

The  Committee  on  Tidal  Hydraulics  conducted  an  aiudysis  of  the  inlet 
during  the  period  March-October  1993,  and  prq)ared  this  rqwrt.  District 
liaison  was  provided  by  Messrs.  Jeffrey  A.  Gebert  and  Keith  D.  Watson. 

The  Committee  on  Tidal  Hydraulics  expresses  its  deq>  gratitude  to 
Messrs.  Gdxrt  and  Widson  for  excellent  briefings,  exhaustive  efforts  in 
locating  and  furnishing  data,  and  personal  msights  into  the  processes  at  Indian 
River  Inlet.  The  Committee  also  thanks  Messrs.  Trimbak  M.  Parchure,  Jeff 
Lillycrop,  and  Allen  M.  Teeter  of  the  U.S.  Anny  Engineer  Waterways 
Experiinem  Station  and  Donald  Raney  of  the  University  of  Alabama  for  their 
prmipt  and  helpful  response  to  questions. 

Mr.  Frank  A.  Hemnaim,  Jr.,  is  Chairman  of  the  Committee  on  Tidal 
Hydraulics,  and  Mr.  Samuel  Powell  is  Headquarters,  U.S.  Army  Corps  of 
E^ineers,  Liaison. 
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Conversion  Factors, 
Non-SI  to  SI  Units  of 
Measurement 
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Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  units 
as  follows: 


MuMply 

■v 

To  Obtain  | 

cubic  fact 

0.02831685 

cubic  maters  | 

cubic  yards 

0.7645S49 

cubic  maters  | 

fact 

0.3048 

meters  | 

square  fact 

0.09290304 

square  maters  | 

square  miles 

2.589998 

square  meters  | 

yards 

0.9144 

meters  | 
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1  Introduction 


Background 

1.  The  U.S.  Army  Engineer  District  (USAED),  Philadelphia,  is  responsi¬ 
ble  for  a  U.S.  Army  Corps  of  Engineers  (USAGE)  navigation  project  at 
Indian  River  Inlet,  Delaware  (Figures  1  and  2.)  Indian  River  bdet  has  experi¬ 
enced  progressive  scour  since  about  1940,  and  the  scour  has  accelerated  since 
the  mid  1970s.  The  Philadelphia  District  has  conducted  studies  of  the  inlet 
scour,  including  a  U.S.  Army  Engineer  Waterways  Experiment  Station  (WES) 
Coastal  Engineering  Research  Center  (CERC)  numerical  model  investigation 
and  an  architect/eqgineer  evaluation  of  potential  remedial  measures. 

2.  In  June  1992,  the  Philadelphia  District  asked  the  USAGE  Committee  on 
Tidal  Hydraulics  to  review  the  Indian  River  Inlet  situation  and  address  several 
questions  posed  by  the  District. 


Purpose 

3.  The  purpose  of  this  rqmrt  is  to  answer  the  following  questions  posed 
by  the  District; 

a.  What  factors  caused  the  accelerated  scour  to  begin  during  the  mid- 
1970s? 

(1)  Does  sand  bypassing  affect  the  scour? 

(2)  Does  the  tidal  prism-inlet  area  theory  ^ly? 

b.  What  course  of  action  is  recommended? 

(1)  Is  monitoring  alone  sufficient  at  present,  or  should  planning- 
design-budgeting  for  channel  stabilization  begin? 

(2)  What  is  the  threshold  for  beginning  stabilization? 

(3)  Are  the  architect/engineer  stabilization  study  findings 
appropriate? 
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c. 


(4)  Is  local  cross-section  control  sufficient? 

(5)  Should  pr«c».  flow,  b.  ««p,ed  and  Ak  twd  amored? 

da^  a„ 


(1)  S'^-'n-ore  hydrographic  survey,  bo  scheduled? 

(2)  What  models  are  needed  to  evaluate  the  „  ki 

stabilization  features?  ^  Problem  and  design 
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igure  2.  Indian  River  Inlet  aerial  photograph  (Courtesy  of  Philadelphia  District) 


2  The  Indian  River  Estuary 
Scour  Problem 


4.  The  information  that  follows  was  taken  from  the  stunmaiy  by  Gebert, 
Watson,  and  Rambo  (1992)  and  from  presentations  made  by  the  Philadelphia 
District  at  the  June  1992  meeting  of  tlM  Committee  on  Tid^  Hydraulics. 

5.  Indian  River  Inlet  is  typical  of  many  East  Coast  barrier  island  inlets. 
Behind  the  inlet  lies  Indian  River  Bay  (about  14  square  miles'  in  surface 
area),  which  is  connected  to  Rehoboth  Bay  (about  the  same  size.)  Figures  1-3 
show  the  bay  and  inlet  system  at  present.  Prior  to  1938,  the  inlet  was 
ephemeral,  opening  in  various  locations  after  a  storm,  then  closing  again  as  it 
filled  with  sediment.  Inlets  dredged  by  the  state  also  closed  quickly.  Between 
1938  and  1940,  the  Corps  of  Engineers  built  parallel  jetties  to  create  a  stable 
S(X)-ft-wide  inlet  that  provided  a  navigation  pass  for  recreational  boats. 

Almost  as  soon  as  the  Federal  project  was  completed,  erosion  of  the  unpro¬ 
tected  interior  inlet  shoreline  began. 

6.  In  1941,  dq)ths  throughout  inlet  were  less  than  20  ft  below  mean 
low  water  (mlw).  The  inlet  scoured  both  its  bed  and  uiq)rotected  banks, 
increasing  the  average  depth  by  about  0.5  ft  per  year.  By  1974  some  holes 
were  deeper  than  40  ft.  From  about  1975  on,  the  scour  rate  increased  to 
about  1  ft  per  year,  and  in  1991  nearly  all  of  the  inlet  was  deeper  than  40  ft 
with  some  holes  over  100  ft  deep.  Scour  of  20  to  80  ft  had  occurred  since 
1941.  One  exception  to  the  overall  scour  is  a  high  spot  700  ft  seaward  of  the 
existing  bridge,  which  is  now  remaining  at  a  constant  depth  of  about  35  ft  or 
shoaling  slightly. 

7.  About  250  ft  of  the  seaward  end  of  the  north  jetty  (about  5  ft  per  year) 
has  been  lost  to  progressive  storm  damage.  During  the  1980s  the  inlet  shore¬ 
line  landward  of  the  jetties  was  stabilized  with  revetments.  Diver  inspections 
of  the  jetty  toes  have  revealed  no  damage  from  channel  scour.  The  north  jetty 
has  been  steadily  losing  stones  from  the  end  since  the  1950s.  Scour  may  ^ve 


'  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI  is  presented  on  page  vi. 
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contributed  somewhat  to  that  loss,  but  the  primary  cause  has  been  storm 
damage.' 

8.  The  original  highway  bridge,  constructed  on  piles  which  significantly 
reduced  the  inlet  cross  section,  was  replaced  during  the  late  1960s  through 
early  1970s  with  a  new  bridge  laiilt  seaward  of  the  old.  The  old  bridge  piers 
were  cut  off,  rather  than  removed  entirely,  and  acoustic  surveys  of  the  inlet 
(Lillycrop  et  al.  1991)  show  evidence  of  the  pile  stubs  and  decking  debris  at 
the  location  of  the  old  bridge. 

9.  Scour  around  the  new  bridge  piers  (two,  located  at  about  the  one-third 
points  across  the  channel)  led  to  1989  placement  of  a  riprap  blanket  3  to  6  ft 
thick  and  300  ft  wide  from  bank  to  ba^  around  the  piers.  Design  velocity 
for  the  ripn^  was  IS  ^s,  and  it  does  not  seem  to  have  sustained  any 
significant  damage  since  installation. 

10.  In  19S7  nourishment  of  the  north  (downdrift)  beach  began,  using 
sediments  from  the  fringes  of  Rehoboth  Bay.  Beginning  in  1972,  sediment 
was  mined  from  the  flood  tide  bar  of  the  iidet  to  nourish  the  beaches.  Mini¬ 
mum  depths  over  the  flood  shoal  were  nearly  zero  before  mining  began.  The 
estimated  flood  shoal  accumulated  sediment  rate  was  about  60,000  cu  yd  per 
year  (USAED,  Ptiiladelphia  1984). 

11.  In  1990,  beach  sand  bypassing  was  begun,  passing  on  average  about 
100,000  cu  yd  per  year  from  Ae  south  beach  to  the  north.  The  estimated 
gross  littoral  transport  is  300,(XX)  to  S(X),000  cu  yd  per  year,  and  the  net 
transport  is  100,000  to  150,000  cu  yd  per  year  to  the  north  (USAED, 
niiladelphia  1984).  (These  are  the  best  estimates  available,  based  on  calcula¬ 
tions  made  using  Wave  Information  Study  (WIS)  data  from  Station  65  north  of 
the  inlet.  Confidence  in  the  results  may  be  reduced  somewhat  by  the  fact  that 
conq)utations  using  the  nearest  WIS  station  (Station  66)  show  a  net  littoral 
transport  to  the  south,  a  result  contradicted  by  all  observations.)  Prior  to 
sand  bypassing,  the  ebb  shoal  accumulated  sediment  at  a  rate  of  about 
100,000  cu  yd  per  year.  Data  since  bypassing  began  are  insufficient  to  tell  if 
that  accumulation  has  continued  (Gebert  et  al.  1992). 

12.  A  chronological  history  of  significant  inlet  events  is  given  below. 


1938 

14-  by  200-ft  Federal  channel  completed  | 

1939 

Jetties  completed 

1940 

Bridge  opened  to  replace  earlier  structure 

1941 

Sheet  pile  bulkhead  constructed  on  inner  shoreline  | 

1947 

Stone  fill  added  to  support  pile  bulkhead  1 

'  Personal  Communication,  30  August  1993,  Jeff  Gebert,  U.S.  Army  Engineer  District, 
Philadelphia,  Philadelphia,  PA. 
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1957 


Outer  MCtion  of  jotties  rapairad 


1967 


Artificial  nourishmant  of  north  beach  begun 


1961 


48,000  cu  yd  mined  from  flood  shoal 


1963 


Sheet  pile  bulkheads  repaired  and  extertdad  westward 


1965 


New  bridge  opertad 


1972 


Removal  of  old  bridge  begun 


1973 


Steal  bulkhead  walls  replaced  by  rubble  mound 


1975 

1976 
1978 
1984 

1989 

1990 


774,000  cu  yd  mined  from  flood  shoal _ 

143,000  cu  yd  mined  from  flood  shoal 
Removal  of  old  bridge  completed 
700,000  cu  yd  mined  from  flood  shoal 
468,000  cu  yd  mined  from  flood  shoal 
Stone  blanket  placed  at  bridge 

1 75,000  cu  yd  mined  from  flood  shoal _ 

Sand  bypassir«g  begun,  100,000  cu  yd  per  year  to  north 
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Inlet  Processes 


Sediments 

13.  Geotechnical  investigations  have  shown  that  the  barrier  island  medium 
to  fine  sands  extend  down  to  about  -35  ft  NGVD,  varying  ftom  -23  ft  at  the 
jetty  tips  to  -40  ft  at  the  south  abutment  of  the  highway  bridge.  Representa¬ 
tive  grain  sizes  {Djj  or  35  percent  finer  by  weigin)  are  0.35  mm  on  the  flood 
shoal,  0.28  mm  on  the  ebb  shoal,  and  0.30  to  0.42  mm  on  the  beaches. 

Below  the  sands  are  lagoonal  clays  and  silts  to  a  depth  of  about  -50  to  -100  ft 
NGVD,  below  which  consolidate  Pleistocene  sands  and  gravels  occur.  Deep 
cores  in  the  area  are  sparse,  so  the  Pleistocene  sediments  may  be  closer  to  the 
surface  in  some  locations. 

14.  The  clays  and  silts  layer  is  soft  at  the  interface  with  the  surface  sands, 
with  a  split  spoon  sandier  requiring  only  two  to  six  blows  per  foot.  Some¬ 
what  lower  in  the  layer  the  sediments  are  slightly  stiffer,  requiring  four  to  six 
blows  per  foot.  KnA  and  John  (1976)  confirm  this  description,  reporting  that 
the  iqqier  portion  (about  half)  of  the  layer  consists  of  soft  lagoonal  sediments, 
and  die  lower  portion  is  made  up  of  somewhat  stiffer  marsh  dqiosits.  Cores 
show  some  thin  layers  of  peat  (USAED,  I%iladelphia  1984). 


Hydrography 

15.  Cross-sectional  dqith  profiles  have  been  regularly  measured  by  the 
Philadelphia  District.  Hy^ographic  data  used  here  are  at  inlet  ranges  spaced 
500  ft  apart  from  1,(X)0  ft  east  (Range  10+00  E)  of  the  bridge  to  1,000  ft 
west  (R^e  10+00  W)  of  the  bridge.  Aj^ndix  A  contains  a  tabular  presen¬ 
tation  of  the  profile  data.  Figure  3  shows  the  range  locations.  Plates  1-5 
show  cross-sectional  profiles  of  the  five  ranges  as  they  have  varied  since 
1941.  Only  Range  0+00  (Plate  3)  exhibited  a  shallower  maximum  depth  in 
1992  than  in  1988,  because  of  the  ripnq>  blanket  placed  at  the  bridge  in  1989. 

16.  Figure  3  shows  the  inlet  dq)ths  for  June  1991  in  contour  form  along 
with  the  range  numbering  system  uUd  in  hydrographic  surveys.  The  inlet 
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exhibits  a  sinuous,  basin-ridge  topography  with  the  following  major  features, 
described  from  the  ocean  landwa^. 

a.  A  basin  (maximum  dq>ths  greater  than  70  ft  below  NGVD)  lies 
just  eastward  of  a  line  connecting  the  ends  of  the  jetties.  The 
centerline  of  this  basin  is  located  in  the  northern  half  of  the  cross 
section. 

b.  A  ridge  (maximum  dq>th  35  ft  below  NGVD)  is  centered  about 
Range  6+8S  E.  The  ridge  saddle  is  toward  the  north  side  of  the 
section. 

c.  Dq)ths  increase  relatively  slowly  westward  of  the  ridge,  with  a 
maximum  depth  of  80  ft  in  a  basin  at  about  Range  1  -l-SS  E.  The 
centerline  of  the  basin  is  slightly  south  of  the  center  of  the  section. 

d.  A  ridge  (maximum  depth  55  ft  below  NGVD)  is  located  at  about 
Range  1+55  W,  corresponding  closely  to  the  location  of  the  old 
bridge. 

e.  A  basin  with  maximum  depths  greater  than  105  ft  below  NGVD 
lies  at  about  Range  4+(X)  W,  south  of  the  centerline  of  the  section. 

/.  Maximum  dq>ths  then  decrease  in  the  westward  direction,  shoaling 
to  less  than  30  ft  below  NGVD  at  about  Range  1 1  +(X)  W. 

17.  Table  1  shows  the  variation  of  the  unweighted  average  inlet  cross- 
sectional  area  from  1941  to  1992  for  the  five  ranges  shown  in  Figure  3.  (The 
basic  data  are  shown  in  Appendix  A,  Tables  A1-A5.)  Range  0+00  data  were 
excluded  from  the  average  for  1992  because  of  the  ripr^  placement.  The 
steady  increase  in  depth  is  obvious,  but  the  data  suggest  that  the  rate  of 
increase  is  slowing. 

18.  Figures  4  and  5  show  the  rate  of  deepening  with  plots  of  overall 
average  inlet  dq)th  and  maximum  inlet  depth  from  1941  to  1992  for  the 
2,000-ft-long  inlet.  Figure  4  shows  that  the  average  inlet  depth  increased 
slowly  (about  0.2  ft  per  year)  from  1941  to  1966,  then  began  to  increase  more 
rigidly  (about  0.8  ft  per  year).  The  rate  of  increase  seems  to  have  slowed 
considerably  in  the  1988-1992  period,  but  the  data  are  too  sparse  to  draw  a 
conclusion  at  this  point.  The  second  curve  on  Figure  4  shows  the  average 
dq>th  of  the  central,  300-ft-wide  portion  of  the  inlet  from  Range  10+00  E  to 
Ri^e  10+(X)  W,  which  excludes  the  protected  flanks  of  the  inlet. 

19.  Figure  5  plots  the  maximum  depth  at  Range  5+00  W  (the  deepest 
point  throughout  most  of  the  period)  and  the  rate  of  depth  increase  at 

Range  5+00  W  over  the  period  of  interest.  Figure  5  shows  a  spurt  in  the  rate 
of  increase  in  maximum  dq>th  after  1974.  Unlike  the  average  depth, 
maximum  dqtth  shows  no  signs  of  leveling  off  after  1988.  It  i^pears  from 
Figures  4  and  5  and  Plates  1-5  that  the  surface  sand  layer  (bottom 
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Irable  1 

1 

lAvmrage  Inlat  Cross-Sactional  Araa,  Range  10  +  00E  to  1 

[Range  10  +  00W 

Area 

Rate  of  Change,  Average  Cross-Seetianai  Area 

Yaw 

sqft 

sqfl/year 

1941 

5,150 

- 

1950 

8,300 

350 

1960 

9,570 

119 

1966 

11,200 

308 

1974 

13,200 

240 

1978 

16,900 

837 

1988 

22,000 

541 

|l992' 

22,900 

200 

I'  Excludes  Range  0+00. 

1 

elevation  -35  ft  NGVD)  began  to  disappear  during  the  late  1960s  and  was 
mostly  gone  by  1980. 

20.  Plate  6  shows  the  change  in  average  dqjth  over  time  for  each  of  the 
five  ranges.  All  five  ranges  exhibit  similar  patterns,  with  the  exception  that 
Range  S+00  W  has  increased  in  depth  much  more  rapidly  than  the  others 
since  1978.  On  average,  the  curves  for  the  other  ranges  suggest  a  decrease  in 
the  rate  of  deepening  after  1978. 

21.  Plate  7  illustrates  longitudinal  profiles  of  average  depth  for  five  of  the 
years  of  record.  It  shows  the  progressive  movement  of  deepening  from  the 
seaward  end  upstream,  culminating  in  the  r^id  erosion  of  Range  S+OO  W 
from  1974  through  1992.  That  sequence  is  consistent  with  scour  of  the  sand 
bed  beginning  at  the  seaward  end  and  progressing  inland  as  the  supply  of  sand 
for  a  hypothesized  net  landward  transport  is  sharply  decreased.  (P^- 
graphs  44-47  show  the  net  landward  hypothesis  to  be  reasonable.)  Placement 
of  the  ripriq)  blanket  at  the  bridge  (Range  0+(X))  in  1989  may  have  con¬ 
tributed  to  the  1988-1992  erosion  at  5+(X)  W  by  partially  blocking  the  sand 
supply  that  would  have  been  contributed  from  the  offshore  area. 

22.  Plate  8  is  a  time-history  plot  of  stations  across  the  channel  at 
Range  5+00  W.  The  stations  indicate  distance  from  the  north  bank  in  feet. 

23.  The  rate  of  inlet  depth  changes  is  displayed  in  thr^  ways  in  Table  2, 
as  a  sediment  loss  rate  in  cubic  yards  per  year,  as  an  average  depth  increase 
rate  in  feet  per  year,  and  as  a  maximum  inlet  depth  increase  in  feet  per  year. 
All  were  computed  from  Range  10+(X)  E  to  Ra^e  10+(X)  W  using  the  range 
depths  in  Table  1,  with  1992  Range  O-r-OO  data  excluded  because  of  the  riprap 
placement. 
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Figure  4.  Average  inlet  depth,  1941-1992 


Figure  5.  Maximum  inlet  depth  changes.  Range  5-t-OOW 
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24.  The  interior  shorelines  of  the  inlet,  landward  of  the  bridge, 
e]q)erienced  severe  erosion  almost  immediately  after  completion  of  the  jetties. 
Keegan  (1967)  reported  that  the  width  at  old  Range  40+00  (sq^roximately 
new  Range  2S+(X)  W)  increased  ftxnn  720  ft  in  1939  to  1,270  ft  in  1948, 
while  the  average  depth  decreased  from  -9  ft  to  about  -8  ft  msl.  In  1989  that 
section  was  more  thsm  2,000  ft  wide  and  had  an  average  dq>th  of  about  20  ft, 
with  dq)ths  farther  inland  very  shallow,  generally  less  than  4  ft. 


Hydrodynamic  Behavior 

25.  Table  3  shows  some  represoitative  tide  ranges  for  the  Indian  River 
Inlet  and  Bay  system.  Station  locations  are  shown  in  Figure  1 . 

26.  Figures  6  and  7  show  tide  measurements  taken  by  WES  in  1988  (Lil- 
lycrop  et  al.  1991),  which  show  the  phase  shifts  and  amplitude  changes  from 
die  ocean  tide  through  the  inlet,  Indian  River  Bay,  and  Rehoboth  Bay.  Figure 
8  shows  that  the  inlet  cross-sectional  avenge  velocity  (at  about  Range  7+00 
E)  leads  the  inlet  tide  by  about  1-1/2  hr,  demonstrating  a  wave  character  inter¬ 
mediate  between  progressive  wave  behavior  and  standing  wave  behavior. 

27.  Raney,  Doughty,  and  Livings  (1990)  performed  a  two-dimensional 
numerical  model  study  of  the  ocean-inlet-bay  system.  They  verified  the  WES 
Implicit  Flooding  Model  (WIFM)  model  to  1988  tide  and  velocity  data,  then 
ran  a  series  of  tests  with  geometry  from  1941  to  1988  plus  several  geonKtries 
with  inlet  size  larger  than  that  of  1988.  The  peak  tidal  discharge  for  the  1941 
geometry  was  1S,0(X)  cfs  at  an  inlet  velocity  of  3  fps,  and  the  peak  1988 
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Figure  6.  Ocean  tide  (adapted  from  Liilycrop  et  al.  (1991) 
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Figure  7.  Measured  interior  tides  (adapted  from  Ullycrop  et  al.  (1991)) 


Figure  8.  Average  ocean  tide  and  inlet  velocity,  1988  (adapted  from 
Ullycrop  et  al.  (1991)) 
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TIDE  LEVEL.  FT  NGVD 


geomeay  disdiarge  was  over  100.000  cfs  at  about  7  They  tested  inlet 
sizes  larger  than  existing  1988  conditions,  including  cross-sectional  areas  of 
about  21.000  to  28.000  sq  ft.  and  found  that  while  tidal  prism  increased 
monotonically.  inld  velocities  peaked  at  about  the  1988  inlet  size  and  then 
declined  with  each  increase  in  inlet  area. 

28.  Keulegan  (1%7)  calculated  a  conqxMite  roughness  of  0.046  for  the 
inlet  for  1948  conditions,  and  Raney.  DougiUy.  and  Livings  (1990)  devel(^>ed 
a  range  of  cell  roughnesses  of  0.03S  to  0.04S  in  verifying  the  WII^  model. 
The  coefficients  ate  not  directly  comparable  because  of  the  difference  in  one- 
and  two-dimensional  schematizations.  but  do  indicate  the  most  probable  range 
of  roughness  coefficient. 


Stable  Inlet  Theory 

29.  The  Jarrett-O’Brien  (Jarrett  1976)  tidal  prism-inlet  area  relationship  is 
an  order  of  magnitude  ^roximation  for  inlm  in  which  the  sediment-scouring 
capacity  of  tidal  flows  is  roughly  in  balance  with  the  suf^ly  of  sand  into  the 
inlet  through  longshore  transport,  wave  suspension,  and  tidal  flows.  (This 
balance  is  sometimes  e}q)ressed  by  the  ratio  of  tidal  prism  to  gross  littoral 
transport.  For  exanq>le,  see  Bruun  (1977).)  This  relationship  cannot  be  used 
to  pi^ct  an  equilibrium  size  for  Indian  River  Inlet  because  of  the  original 
serious  imbalance  in  sediment-scouring  power  versus  supply,  and  due  to  the 
sand  bypassing  q)eration  which  now  makes  its  controlling  conditions  deviate 
further  from  that  of  otherwise  similar  inlets.  Nevertheless,  it  is  instructive  to 
examine  what  the  relationship  says  about  Indian  River  Inlet. 

30.  Lillycn^  et  al.  (1991)  performed  a  hydrodynamic  analysis  using  an 
adaptation  of  Keegan’s  (1%7)  one-dimensional  c^culation  for  inlet-bay  tidal 
hydraulics.  (The  latter  was  an  effort,  sponsored  by  the  Committee  on  Tidal 
Hydraulics,  that  used  Indian  River  Inlet  as  an  example.)  The  results  showed 
that  the  inlet  could  reach  a  stable  size  of  about  40,000  to  80,(X)0  sq  ft  when 
the  Keulegan-King  hydrodynamics  curve  imersected  the  Jarr^-O’Brien  stabil¬ 
ity  curve.  The  latter  result  can  also  be  obtained  by  the  short  cut  of  assuming 
that  tide  ranges  inside  Indian  River  and  Rehoboth  Bays  are  both  equal  to  the 
ocean  tide  spring  range,  then  using  that  range  to  conqnite  a  tidal  prism  that  is 
used  in  the  Jarrett-O’Brien  equation. 

31.  A  somewhat  more  realistic  maximum  size  can  be  estimated  by  cal¬ 
culating  n  tidal  prism  for  ranges  of  S.3  ft  in  Indian  River  Bay  and  2.65  ft  in 
Rehoboth  Bay,  based  on  the  relationship  between  the  two  bays’  tides  shown  in 
Table  3.  That  calculation  yields  a  tidal  prism  of  3.2  x  10^  cu  ft  and  a  cross- 
sectional  area  of  about  62,0(X)  sq  ft.  The  inlet  would  have  to  have  an  average 
dq>th  of  124  ft  for  that  cross-sectional  area. 

32.  Another  limitation  of  the  Jarrett-O’Brien  equation  for  this  purpose  is 
that  the  scatter  in  the  raw  data  is  large.  Using  the  9S-percent  confidence 
limits  on  that  equation  and  a  prism  of  3.2  x  10^  cu  ft  yields  a  range  of 
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proi>able  cross-sections  of  about  25,000  to  140,000  sq  ft,  or  a  mean  depth  of 
SO  to  280  ft.  a  range  which  is  hardly  amenable  to  arriving  at  a  decision  on 
protective  works. 

33.  O’Brien  and  Dean  (1972)  presented  a  useful  metlKxl  of  plotting  maxi¬ 
mum  inlet  velocity  versus  inlet  cross-sectional  area  for  a  given  inlet,  then 
noting  thitt  certain  parts  of  the  curve  are  inherently  stable  or  unstable.  Such  a 
curve  (see  paragraph  39)  was  construaed  for  Indian  River  Inl^  using  the  data 
from  the  WIFM  n^el  plus  results  from  an  application  of  Keulegan’s  (1967) 
tidal  repletion  analysis. 

34.  A  direct  extension  of  Keulegan’s  computations  for  the  present  Indian 
River  Inlet  was  not  possible.  Keegan’s  original  analysis  used  an  inlet  length 
of  5,500  ft  and  an  equivalent  prismatic  inlet  area  of  9,660  sq  ft.  Erosion  of 
the  interior  shoreline  has  now  made  that  a  very  poor  ^roximation  to  the 
system  since  the  defined  inlet  is  now  only  about  2,000  ft  long,  and  the  rest  of 
the  former  inlet  is  difficult  to  distinguish  ftom  the  bay  pn^r.  Using  2,000  ft 
for  inlet  length,  and  a  weighted  average  area,  a  Manning’s  n  of  0.08  was 
found  to  be  necessary  to  rq>roduce  the  observed  1988  tides  of  Figure  7,  which 
is  clearly  not  a  reasonable  value. 

35.  Although  the  interior  portion  of  the  former  inlet  is  much  wider  now,  it 
must  still  exercise  some  control  over  tides  in  the  system  through  the  extremely 
shallow  flood  tide  shoal;  therefore,  that  effea  was  included  by  treating  it  as  an 
entrance/exit  loss.  Keulegan's  (1967)  nondimensional  head  loss  term  H  was 
of  the  form 


aH  =  ^  *  m  (1) 

R 


where 

f  =  inlet  Darcy  friction  factor 
L  =  inlet  length 
R  =  inlet  hydraulic  radius 
m  =  entrance/exit  loss  coefficieiU,  assumed  =  1 

Another  loss  term  was  added  to  the  equation  to  account  for  head  loss  across 
the  flood  tide  shoal,  using 

aH  =  *  aH'  (2) 

R 


where 
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f  =  Darcy  friction  factor  for  the  shoal 

L'  =  lengtii  of  shoal  zone 

R'  =  hydr»ilic  radius  of  the  shoal 

36.  The  primed  term  for  the  shoal  head  loss  is  similar  to  that  of  the  inlet 
friction  term,  but  it  was  treated  as  a  lumped  head  loss,  using  the  form  to  com¬ 
pute  a  range  of  reasonable  values.  Darcy  coefficients  rapging  from  0.03  to 
0.07,  an  L'  of  2,000  to  3,500  ft,  and  an  i?'  of  10  to  12  ft  yields  a  aH’  value 
of  6  to  25.  By  trial,  tide  ranges  were  confuted  as  descriM  below  to  find  a 
value  of  aH'  that  gave  an  Indian  River  Bay  tide  range  equal  to  that  measured 
in  1988.  This  method  (ignormg  Rehoboth  Bay)  will  give  a  tidal  prism  that  is 
slightly  too  low,  but  close  enough  for  the  present  purpose. 

37.  For  the  tidal  prism  calculations,  2,000  ft  was  used  for  the  inlet  length 
and  a  uniform  500  ft  for  the  width.  For  1988  and  previous  years,  the 
representative  inlet  cross-sectional  area  was  based  on  a  weighted  average  of 
the  five  cross  sections,  neglecting  those  portions  of  5-1-00  W  and  10-1-00  W 
that  lay  outside  the  500-ft-wide,  deep-water  portion.  (These  confutations  do 
not  require  that  the  inlet  be  of  constant  cross  section,  but  only  that  the  cross 
section  be  uniform  enough  to  permit  a  calculation  of  average  flow  to  reason¬ 
ably  rq)resent  average  conditions.  The  more  nonuniform  the  inlet,  the  less 
accurate  the  results  become.  As  discussed  later,  an  inf  rovement  to  these 
simple  calculations  requires  detailed  modeling.)  The  bay  area  was  4.2  x  10* 
sq  ft.  A  Manning’s  n  value  of  0.035  was  used  for  the  inlet.  A  lumped  flood 
shoal  loss  term  aH'  =  14  successfully  reproduced  1988  observed  tides,  and 
that  value  was  used  throughout. 

38.  The  inlet  areas  and  resulting  tidal  prisms  from  the  Keulegan  method 
are  shown  in  Table  4,  along  with  a  calculated  maximum  inlet  velocity  based 
on  an  assumed  sinusoidally  varying  velocity  time-history.  The  sinusoidal 
assunf  tion  is  a  particularly  rough  one,  and  should  be  used  with  caution. 
Shown  also  are  the  maximum  calculated  velocities  from  the  WIFM  model, 
which  enfloyed  a  larger  tide  range,  5.3  ft,  than  the  4.1-ft  average  spring 
range  employed;  therefore,  the  velocities  were  adjusted  by  the  ratio  of  tide 
ranges.  The  velocities  are  conf  arable,  and  both  show  that  peak  velocities 
occur  for  a  cross-sectional  area  near  the  1992  condition.  For  further  inlet 
enlargement,  the  tidal  prism  continues  to  grow,  but  the  inlet  area  grows  faster, 
so  maximum  velocities  begin  to  fall. 

39.  The  data  of  Table  4  are  presented  in  an  O’Brien-Dean-type  plot  in 
Figure  9.  Also  plotted  are  a  maximum  velocity  form  of  the  Jarrett-O’Brien 
equatkm  and  two  field  measurement  points  from  1941  (Keulegan  1967)  and 
1988  (Lillycrop  et  al.  1991).  The  V^FM  velocities  and  the  1988  field 
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observation  have  been  adjusted  (multiplied  by  the  tide  range  ratio)  to  make 
them  conq)arable  with  the  4. 1-fit  tide  range  employed.  In  Figure  9,  the 
Jarr^-O’Brien  curve  represents  the  line  along  which  all  Atlantic  coast  inlets 
with  two  jetties  tend  to  lie,  whereas  the  other  curves  represent  the  specific 
response  of  Indian  River  Inlet. 

40.  The  most  obvious  fact  demonstrated  by  Figure  9  is  that  the  Jarrett- 
O’Brien  curve  means  that  stable  inlet  velocities  for  the  Atlantic  coast  are 
generally  between  3  and  3.S  fjps,  far  less  than  the  present  normal  spring  range 
Indian  River  Inlet  velocities  of  more  than  4.5  Indian  River  Inlet 
velocities  have  not  been  near  that  equilibrium  range  since  the  1940s.  The 
inlet  velocity  was  below  the  Jarrett-O’Brien  line  in  1948,  but  was  obviously 
unstable.  Tlie  velocities  for  1988  geometry  (about  20,000-sq-ft  area)  are,  as 
observed  in  Table  4,  patently  near  the  maximum  expected,  and  may  even 
decrease  with  further  erosion.  The  computed  Indian  River  Inlet  velocities  are 
highly  a{q>roximate  for  geometries  and  roughness  far  from  the  verification 
condition,  but  they  do  illustrate  the  basic  situation:  for  inlet  sizes  larger  than 
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Figure  9.  Maximum  inlet  velocity  as  a  function  of  cross  section 

those  at  present,  the  velocities  may  be  flat  or  even  decreasing,  but  they  are 
still  far  above  the  stability  level. 

41.  The  sht^  of  the  Keulegan  method  and  WIFM  curves  illustrates  an 
inqmrtant  aspect  of  the  O’Brien-Dean  curve.  Points  along  the  rising  side 
(areas  less  than  about  20,000  sq  ft)  tend  to  be  unstable  because  changes  in 
area  tend  to  be  anq>lified;  a  decrease  in  area  reduces  velocity,  which  further 
reduces  sediment-scouring  capacity  and  leads  to  still  more  inlet  area  reduction, 
and  inlet  area  increase  acconq>lishes  the  reverse.  For  points  to  the  right  (fall¬ 
ing  side)  of  20,(XX)  sq  ft,  changes  in  area  cause  compensating  changes  in 
velocity,  so  either  erosion  or  infilling  is  opposed.  These  are,  however,  just 
tendencies  toward  stability  or  instability,  which  ultimately  depends  on  both  the 
hydrodynamic  response  shown  in  Figure  9  plus  the  supply  of  sand  to  the  inlet. 

42.  Figure  9  could  be  used  to  support  an  equilibrium  inlet  size  of  about 
34,(XX)  sq  ft,  the  point  at  which  the  Keulegan  curve  intersects  the  Jarrett- 
O’Brien  curve;  however,  that  cannot  be  made  a  conclusion  at  this  point. 

First,  the  Jarrett-O’Brien  curve  is  based  on  average  sand  supply  conditions  for 
the  Atlantic  coast,  and  it  appears  the  Indian  River  Inlet  supply  is  quite  small. 
Second,  as  noted  above,  the  data  scatter  in  the  Jarrett-O’Brien  curve  is  large. 
Third,  the  Keulegan  method  velocity  curve  is  highly  approximate  and  does  not 
adequately  include  the  Rehoboth  Bay  response.  Figure  9  definitely  indicates 
that  a  stable  sand  bed  is  not  expected  under  existing  conditions,  even  without 
sand  bypassing  in  operation. 
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Sand  Transport 


43.  As  a  first  step  to  identifying  the  stability  of  the  inlet,  a  set  of  sand 
transport  calculations  was  made  to  identify  the  size  of  the  sand  supply  deficit. 
The  calculations  were  performed  as  follows: 

a.  The  inlet  velocities  measured  in  1988  were  smoothed  such  that  a 
single,  repetitive  12.42'hr  tidal  cycle  of  hourly  velocities  was 
obtained,  while  keq)ing  the  integral  over  time  of  the  velocity  and 
velocity  squared  consistent  with  the  raw  values.  This  provided  a 
velocity  time-histoiy  that  preserved  the  distortion  of  inlet  velocities 
by  geometric  effects. 

b.  The  velocity  history  from  a  was  scaled  down  to  average  inlet 
velocities  by  multiplying  the  individual  values  by  three  coefficients 
as  shown  below. 


C,  *Cj*C3*V,^ 


(3) 


where 

C,  -  the  ratio  of  cross-sectional  area  where  the  measurements  were 
taken  (15,018  sq  ft)  to  an  average  cross-sectional  area  Ay^  for 
the  year  (10  and  location  (L)  within  the  inlet  needed 

Cj  =  the  ratio  of  the  ocean  tide  range  for  which  velocities  were 
desired  (e.g.,  4.1  ft)  to  that  of  the  1988  measurement  period 
(5.3  ft) 

Cj  =  the  Indian  River  Bay  tide  modulation  factor 


C3 


(4) 


R  =  tide  range  in  Indian  River  Bay  for  an  ocean  tide  range  of  4.1  ft 
(mean  spring  tide),  and  the  subscript  1988  indicates  the  year. 

The  velocities  were  corrected  with  coefficient  Q  for  spring,  mean, 
and  neq)  tide  ranges  as  listed  for  the  tide  tables  in  Table  3.  (Those 
representative  tide  ranges  were  estimated  since  the  exact  magnitudes 
are  not  critical  to  the  purposes  of  the  calculation.)  The  Q  correction 
was  made  to  account  for  the  change  in  relative  tide  modulation 
between  the  year  1988  and  both  earlier  and  future  cross-sectional 
areas,  as  calculated  by  the  Keulegan  method  and  presented  in  Table  4. 
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c.  The  adjusted  velocity  time-history  was  used  to  compute  plane  bed 
shear  stress  history  by  a  Darcy-type  equation  with  the  friction 
coefficient  derived  from  the  sediment  grain  size  via  Strickler’s 
equation,  and  total  shear  stress  by  Manning’s  equation. 

d.  The  shear  stresses  computed  by  c  were  used  with  sediment 
characteristics  and  flow  information  to  calculate  hourly  rates  of 
sediment  transport  via  the  Ackers- White  total  load  equation 
(Ackers  and  vhiite  1973)  for  one  12.42-hr  tidal  cycle  under  mean, 
average  neap,  and  average  spring  tide  range  conditions. 

e.  The  hourly  sediment  transport  rates  were  integrated  over  a  tidal 
cycle  by  Simpson’s  rule  to  obtain  an  estimate  of  net  and  gross 
sediment  transport  rates  through  the  inlet  for  each  tide  range 
condition. 

/.  The  individual  tide  range  annual  rates  were  linearly  combined  (by 
assuming  that  mean  tide  conditions  predominate  half  the  time,  aid 
typical  ne^  and  spring  range  conditions  each  predominate  one 
quarter  of  the  time)  to  produce  average  annual  net  and  gross  sand 
transport  rates  through  The  inlet. 

44.  The  sand  transport  calculations  used  the  information  above  plus  a 
grain  size  of  0.35  mm  and  Manning’s  n  of  0.035  to  produce  the  said 
transport  results  of  Table  S. 

45.  These  results  must  be  uiderstood  to  be  highly  approximate  (they  do 
not  include  wave  or  storm  effects,  which  will  be  very  significant)  aid  also  to 
be  only  the  sand  transport  capacity.  For  example,  even  though  the  1988  con¬ 
ditions  provide  a  capacity  to  transport  about  5  million  cu  yd  per  year  in  both 
flood  aid  ebb  directions,  with  a  net  flood  direction  transport  of  200,000  cu  yd 
per  year,  the  available  supply  in  the  bed  aid  from  the  entrance  was  much  less, 
so  only  about  60,000  cu  yd  per  year  was  been  transported  into  the  flood  shoal 
area  for  that  period.  (Paragnqih  10  cited  accumulation  of  60,000  cu  yd  per 
year  on  the  flood  shoal,  which  will  be  somewhat  less  than  the  net  transport, 
depending  on  trap  efficiency  of  the  shoal.)  The  results  also  do  not  reflect 
extreme  conditions,  e.g.,  the  largest  tide  considered  is  the  average  spring 
range  of  4.1  ft. 

46.  Inspection  of  the  history  of  depth  changes  along  with  the  thickness  of 
sand  bed  in  the  inlet  led  to  the  conclusion  that  the  sand  bed  within  the  narrow 
portion  of  the  inlet  was  depleted  in  the  mid-1970s.  If  a  supply  from  the  ocean 
of  about  60,000  to  100,000  cu  yd  per  year  is  assumed,  then  Table  5  shows 
why  the  rate  of  bed  loss  continued  at  a  small  rate  (about  15,000  cu  yd  per 
year  from  1941  through  1966,  as  shown  in  Table  2)  until  the  1970s,  when 
sand  transport  began  to  rise,  quickly  stripp'ng  the  bed  of  sand  and  exposing 
the  clay  beneath.  It  is  probable  that,  despite  the  1988  edacity  to  transport  on 
the  order  of  200,(X)0  cu  yd  per  year,  the  actual  rate  was  near  the  60,000  cu  yd 
per  year  observed  to  be  accumulating  on  the  flood  shoal.  That  amount  of 
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supply  would  not  protect  the  bed;  in  fact,  it  would  accelerate  the  scour  of  clay 
by  abrading  the  clay. 
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47.  Table  S  results  suggest  that  for  inlet  sizes  larger  than  about 
30,000  sq  ft,  only  a  relatively  small  quantity  of  sand  would  be  required  to 
maintain  a  bed.  The  inlet  could  well  be  stable  under  the  assumed  conditions. 
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Clay  Erosion 

48.  Clay  erosion  calculations  were  made  to  ex.  nine  the  recent  rates  of 
inlet  enlargement  and  to  characterize  the  bed  change  for  predictions  of  future 
evolution.  Steps  a-c  from  the  sand  transport  calculations  were  again  used, 
except  that  only  the  toud  shear  stress  was  computed,  using  Manning’s  equa¬ 
tion  with  n  again  set  to  0.035.  For  step  d  the  clay  erosion  rates  were 
calculated  by  the  Ariathurai-Partheniades  equation  (Ariathurai,  MacArthur, 
and  Krone  1977), 


E  =  Af(T/T,  -  1) 


(5) 


where 

M  =  erosion  rate  constant,  g/(sq  m  -  min) 
r  =  bed  shear  stress,  N/sq  m 

r,  =  critical  shear  stress  at  which  erosion  begins,  N/sq  m 
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Values  for  hi  and  r,  are  determined  by  flume  tests  and  vary  considerably  fr(»n 
one  cdbesive  sedimeitt  to  another.  Therefore,  fixing  reasonable  values  a  priori 
is  less  straightforward  than  fixing  coefficients  for  sand  transport.  Ranges  of 
reported  values  for  r,  are  given  in  Table  6. 


TaUne 

Rang#  of  Vakias  for  Eroaion  Equation  Coaffidants 

Type  Clays 

»e 

N/sa  m 

Very  soft 

2 

Soft  to  medium  stiff 

7-9 

Stiff 

20 

Values  of  the  erosion  rate  constant  can  range  ff(»n  30  to  300  g/sq  m-min, 
with  the  higher  values  applying  to  stiffer  sediments,  but  without  a  clear 
correlation  to  easily  measu^  parameters  (Ariathurai  and  Arulalandan  1978). 
In  cases  where  sand  is  transported  over  a  clay  bed,  abrasion  by  the  sand 
grains  can  cause  a  somewhat  stiff  clay  to  ero^  as  if  it  were  nnich  softer,  so 
that  the  lower  critical  shear  stress  and  higher  erosion  rate  constants  iq>ply.' 

49.  The  conqtutational  procedure  for  st^s  e  and /was  again  identical  to 
that  described  above  for  sand  transport,  and  clay  erosion  rates  were  calculated 
for  neap,  mean,  and  spring  tides  a^  then  conq}osited  for  a  typical  year  as 
describe  above  for  s^. 

50.  The  sediment  bulk  wet  density  was  d^ermined  to  be  1.70  g/cu  cm 
from  lab  analysis  of  core  sanq>les-S4.9  percent  moisture  content,  67.4  Ib/cu  ft 
dry  density  (USAED,  Philadelphia  1984).  Coefficients  r,.  and  ilf  were  varied 
by  trial  to  generate  erosion  rates  close  to  the  average  and  maximum  rates 
shown  in  Table  2,  resulting  in  values  of  S.S  to  6.6  Msq  m  and  35  to  60 
gm/sq  m-min,  respectively.  A  summary  of  clay  erosion  results  is  given  in 
Table  7. 

51.  The  Table  7  results,  which,  like  the  sand  transport  results,  are  highly 
iqjproxiinate,  suggest  that  the  clay  erosion  will  be  eveiWlly  self-limiting. 

The  clay  will  not  likely  become  more  erodible  with  dq)th,  and  will  likely 
cease  eroding  altogether  at  an  inlet  size  of  about  26,000  to  30,000  sq  ft. 


'  Personal  communication,  2S  May  1993,  Allen  M.  Teeter,  WES,  Vicksburg,  MS. 
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4  Causes  of  Erosion 


52.  From  die  above  analyses,  smne  todative  conclusions  about  the  causes 
of  Indian  River  Inlet  erosion  can  be  drawn;  however,  firm  conclusions  on  the 
quantitative  comribution  of  secondary  factors  are  limited.  Three  events 
occurred  in  the  mid-1970s  roughly  coincideiU  with  an  acceleration  in  bed 
erosion:  the  sand  layer  in  the  inlet  was  stripped  away,  exposing  the  clay 
layer;  the  old  bridge  piles  were  cut  off,  prob^ly  redoing  the  inlet’s  frictional 
resistance  to  flow;  and  the  flood  tide  shoal  was  mined  for  sediment  to  be  used 
in  beach  nourishment.  An  argument  can  be  made  linking  any  or  all  three  of 
these  events  to  the  erosion,  but  the  available  data  do  not  permit  separation  of 
their  relative  contributions. 


Initial  Dredging  and  Jetties 

53.  Once  the  inlet  size  was  established  with  jetties  to  resist  closure,  inl^ 
deq)ening  was  inevitable.  Figure  9  and  Table  5  taken  together  show  that  the 
early  inlet  was  on  the  unstable  portion  of  its  response  curve.  The  jetties 
restricted  the  stqiply  of  sand,  tte  enlarging  tidal  prism  increased  the  transport 
rate  through  the  inlM,  and  scour  commraiced,  puling  the  inlet  \q)  the 
response  curve.  With  width  constrained,  deepening  was  ordained. 


Flow  Convergence  Zone  Scour 

54.  The  deq)  holes  shown  in  Figure  3  are  probably  the  consequence  of 
local  concentration  of  hydraulic  stresses,  perhaps  at  convergence  zones  of 
turbuloit  structures  generated  by  the  geometry /roughness  or  shed  by  obstruc¬ 
tions  such  as  the  new  bridge  piers  and  renmants  of  the  old  bridge. 


Scour  of  Cohesive  Sediment  Beds 

55.  The  rapid  scour  since  the  niid-1970s  is  at  least  in  part  due  to  exposure 
of  the  clay  bed  when  the  sand  layers  were  scoured  away,  and  that  may  in  fact 
be  the  most  significant  cause  of  rapid  erosion.  Clay  eroded  from  the  bed  is 
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immedutely  transported  through  the  inlet  and  will  not  redq)osit  there  in 
significant  quaittities:  therefore,  unlike  sand,  once  clay  is  entrained,  it  is  lost. 
The  nodwinogeneous  nature  of  the  clay  beds  (deq>  cores  are  fairly  ^>arse 
along  the  inlet  centerline)  could  be  a  factor  in  the  development  of  scour  holes 
much  deeper  dian  the  avenge  dqpth,  but  the  flow  convergence  mrationed 
above  is  probably  the  major  contribittor  to  rapid  clay  erosion. 


Old  Bridge  Removal 

56.  Removal  of  the  old  bridge  piers  may  have  reduced  the  inlet’s  effective 
rou^iness  (from  a  Manning’s  n  of  0.046  to  about  0.035)  and  thus  contributed 
to  the  growing  tidal  prism  and  troid  of  increasing  erosion.  However,  the 
evidence  is  not  conclusive,  as  described  above.  The  old  bridge  debris  (pile 
stubs  and  decking)  may  have  localty  decreased  erosion  by  armoring,  causing 
die  ridge  observed  at  Range  1  +55  W  (see  Figure  3)  even  while  contributing 
to  overall  erosion  by  increasing  the  flow  capacity  of  the  inlet. 


Interior  Shoreline  Erosion 

57.  In  1948  the  inlet  was  5,500  ft  long.  Subsequent  interior  shore  erosion 
opened  the  landward  aid  of  the  inlet  to  s^re  the  effective  length  is  now  only 
about  2,000  ft.  That  shoreline  erosion  probably  contributed  to  a  more  nqiid 
increase  in  tidal  prism,  and  thus  more  i^a  erosion. 


Sand  Mining 

58.  Sand  mining  may  have  accelerated  the  erosion  somewhat,  but  the  sand 
transport  data  of  Table  5  suggest  it  was  not  a  critical  factor.  The  calculations 
showed  a  na  landward  sand  transport  rate,  which  suggests  that  sand  deposited 
on  the  shoal  remained  there  and  did  not  reenter  the  inla  under  ordinary  cir¬ 
cumstances.  Increased  wave  penetration  may  have  resulted,  leading  to  more 
interior  shore  aosion,  which  could  have  been  an  indirect  contribution  to  inlet 
erosion. 


Sand  Bypassing 

59.  Sand  bypassing  may  have  increased  the  erosion  rate  somewhat.  Given 
that  the  sand  bypassing  is  about  equal  to  the  na  tidal  transport  in  the  inla  and 
less  than  one  fiftieth  the  gross  transport  cqiacity  within  the  inla,  putting 
100,(X)0  cu  yd  of  sand  per  year  into  the  inla  or  taking  it  out  does  not  seem  to 
be  a  critical  factor  at  this  point. 
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5  Predicting  the  Probabie 
Equilibrium  inlet  Size 


Stable  Sandy  Inlet 

60.  The  Jarrett-O’Brien  curve  will  not  serve  to  predict  an  Indian  River 
Inlet  equilibrium  size.  A  more  appropriate  calculation  is  the  O’Brien-Dean 
stability  curve  in  Figure  9,  which  shows  that  the  inlet  was  unstable  for  sizes 
less  dian  about  20,000  sq  ft  and  tends  toward  stability  for  cross  sections  larger 
dian  that.  Combining  thtt  curve  with  a  velocity  form  of  the  Jarrett-O’Brien 
curve  suggests  that  the  inlet  might  stabilize  a  sand  bed  at  about  34,000  sq  ft,  if 
die  littoral  supply  were  sufficient,  but  it  probably  is  not,  particularly  in  light 
of  the  present  sand  bypassing.  Therefore,  neithn  curve  can  be  used  to 
reliably  predict  the  e^ibrium  size. 


Erosion  Resistance  of  Clay  Sediments 

61 .  The  clay  erosion  calculations  suggest  that  erosion  will  be  self-limiting 
at  an  inlet  size  of  less  than  30,000  sq  ft.  Spots  of  continued  erosion  could 
still  contitBie.  A  cross-sectional  area  of  30,000  sq  ft  will  result  ftom  an  aver¬ 
age  depth  of  60  ft,  which  suggests  maximum  depths  of  twice  that,  or  120  ft. 
However,  the  Pleistocene  sediments  at  about  -100  ft  NGVD  will  probably 
limit  the  maximum  depth. 


Maximum  Scour  Hole  Size 

62.  The  maximum  possible  scour  depth  is  likely  to  be  about  100  ft-the 
depth  of  the  Pleistocene  sands  and  gravels.  At  that  depth  the  sand  transport 
calculations  suggest  a  diminished  flow  transport  cecity  meeting  a 
significandy  less  erodible  material.  Note,  however,  that  this  is  not  a  firm 
conclusion,  since  the  sand  transport  calculations  assume  a  near-uniform 
geometry  and  do  not  consider  flow  convergence  zones  at  all. 
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63.  The  lateral  size  of  the  scour  holes  cannot  be  safely  estimated  by  the 
simple  methods  used  herein,  except  in  the  limit  of  overall  enlargement  of  the 
inlet  to  a  nearly  uniform  dq)th  as  described  in  paragraphs  42-51 . 

64.  The  Cmnmittee  on  Tidal  Hydraulics  notes  that  experience  shows  that 
scour  holes  in  general  seem  to  stabilize  at  dqiths  of  150  ft  or  less;  very  few 
such  holes  ever  go  deeper  than  that.  The  evidence  for  this  limiting  depth  is 
anecdotal  and  undocumented,  but  the  observation  is  supported  by  several 
Cmnmittee  members’  e)q)erience. 
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6  Remedial  Measures 


Architect/Engineer  Recommendations 

65.  The  Philadelphia  District  obtained  a  consulting  firm  to  identify 
possible  stabilization  schemes  (Tetra  Tech  1992)  and  received  the  following 
suggestions  (Committee  comments  in  brackets): 

a.  Rubble  weir  at  Station  7+00  E  (the  shallow  platform). 

[Committee  Comment:  A  single  weir  is  unlikely  to  impede  the  flow  enough  to 
stabilize  or  reduce  the  tidal  prism,  and  it  is  likely  to  accelerate  erosion 
elsewhere.  The  presence  of  ridges  as  described  in  paragn4>h  16  has  not  con¬ 
trolled  the  scour,  and  may  even  be  a  contributing  factor  to  the  localized  holes 
as  mentioned  in  par^raph  56.  Continued  deepening  of  the  adjacent  scour 
holes  could  cause  a  weir  to  collapse.] 

b.  Rubble  weir  at  Station  0  (the  bridge). 

[Committee  Comment:  Unlikely  to  help  for  the  same  reasons  as  given  above. 
The  bridge  protection  riprap  luu  essentially  accomplished  this  already,  though 
the  center  portion  is  lower  than  the  sides.] 

c.  Armor  the  entire  inlet. 

[Committee  Conunent:  It  would  work,  but  the  enormous  cost  makes  it  a  last 
resort.  It  is  probably  overkill,  in  that  some  lesser  solution  will  probably 
suffice.] 

d.  Line  the  jetty  toes  with  erosion  control  mats  (HYDRACOR  mats 
with  froruilike  extensions  on  the  surface). 

[Committee  Comment:  Might  help,  but  at  present  the  jetty  toes  are  not  threat¬ 
ened.  If  they  are,  another  aiproach,  such  as  launching  stone  techniques 
(Headquarters,  USACE  1991),  may  also  be  useful.] 

e.  A  submerged  sill  with  erosion  control  mats  at  station  7+00. 
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[Committee  Comment:  Unlikely  to  belp.  Tbe  comments  on  item  a  also  apply 
to  e.] 


f.  Submerged  groins  at  station  7+00  E. 

[Ctmonittee  Comment:  Unlikely  to  help.  The  comments  on  item  a  also  apply 
to/.] 


g.  A  wood  pile  weir  250  ft  inland  from  the  bridge. 

[Committee  Omunent:  Based  on  the  hydraulic  calculations  here,  adding 
roughness  will  not  stabilize  the  tidal  prism.  As  stated  above,  a  single  weir  at 
ai^  location  is  not  likely  to  help.] 

66.  Some  coiid>ination  of  these  measures,  such  as  weirs  spaced  at  regular 
intervals,  could  be  effective  if  designed  properly.  Extensive  design  studies 
would  be  required  for  any  such  plan  in  ordn  to  test  its  effectiveness  and  lack 
of  harmful  effects.  Spaced  weirs  mi^  be  so  close  together  that  they 
approadied  continuous  armoring. 


Other  Possible  Remedial  Measures 

AdcHtional  inlot 

67.  An  additional  inlet  would  comribute  part  of  the  tidal  prism,  perluqps 
reducing  the  flow  throu^  Indian  River  Inlet  and  lessening  the  erosive  pres¬ 
sure.  A  direct  ocean  connection  to  Rdtoboth  Bay  would  also  contribute  to 
improved  circulation  and  flushing,  possibly  improving  water  quality  in  both 
bays.  (An  extremely  rigorous  design  study  wc^d  be  necessary  to  obtain  the 
potential  benefits.)  However,  cutting  a  new  inlet  would  further  interrupt  the 
littoral  sand  stq^ly  in  a  region  already  experiencing  beach  erosion.  Public 
and  political  q)proval  of  another  inlet  is  unlikely. 


Offshoro  sand  mining  for  inlat  nowishment 

68.  It  might  be  possible  to  protect  the  clay  bottom  from  further  erosion  by 
maintaining  a  carpet  of  sand  movii^  across  the  inlet  bottom.  The  basic  pro¬ 
cess  would  consist  of  occasionally  dredging  suitable  sand  firom  offshore,  plac¬ 
ing  it  in  the  inlet  throat,  and  then  recycling  it  by  regularly  mining  from  the 
flood  shoal  and  reintroducing  it  to  the  inlet.  A  few  hundred  thousand  cubic 
yards  would  cover  the  bottom,  but  the  sand  transport  rates  of  Table  5  suggest 
that  several  million  cubic  yards  of  sand  would  be  required  to  keep  the  inlet 
safely  covered.  Once  the  bed  was  covered,  the  recycling  of  sand  from  the 
flood  shoal  would  probably  have  to  be  at  least  2(X),000  cu  yd  per  year.  While 
this  potential  solution  has  several  merits,  including  potentially  contributing  to 
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beach  oouridunent  by  the  amount  of  sand  expeUed  by  each  d)b  {riiase,  it  is 
unlikdy  to  be  an  ectmomical  solution. 


Sand  wair  in  jattias 

69.  It  has  been  q)eculated  that  the  ridge  observed  in  Figure  3  at  about 
Range  6+85  E  is  caused  by  sand  leaking  through  and  over  the  south  jetty  and 
providing  a  protective  layer  that  inhibits  erosion  (much  like  the  sand  nourish¬ 
ment  described  above).  If  that  is  happening,  then  assisting  the  process  by 
creating  a  path  for  sand  flow  at  that  point  could  reduce  the  rate  of  erosion 
adjacoit  to  the  ridge.  Such  a  path  c^d  be  created  by  lowering  the  jetty  to 
mlw  as  in  weir  jetty  construction,  or  by  sinq>ly  removing  a  few  jetty  stones  to 
make  the  jetty  more  porous.  Consideration  of  such  a  step  will  require  careful 
consideration  of  the  possibility  that  it  will  exacerbate  beadi  erosion  south  of 
the  inlet. 
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7  Conclusions  and 
Recommendations 


70.  The  data  and  analyses  presented  here  suggest  that  overall  inlet  erosion 
may  be  slowing.  Tlie  data  do  not  suggest  that  deep  hole  scour  is  slowing,  but 
it  is  believed  that  the  deep  hole  at  Range  S+00  W  may  not  deepen  further 
since  it  has  encountered  the  Pleistocene  sediments  below  the  clay. 

71 .  The  scour  does  not  yet  present  a  clear  danger  to  the  bridge  and  jetties 
that  justifies  remedial  measures  at  this  time.  Continued  regular  monitoring  is 
clearh  essential,  with  contingency  plans  in  place  to  begin  action  if  remedial 
action  is  needed. 

72.  A  geotechnical  analysis  of  slope  stability  to  predict  at  what  scour  hole 
size  the  jetties  or  bridge  will  be  endangered  is  recommended.  The  analysis 
will  provide  a  criterion  against  which  to  compare  future  scour  hole  growth 
and  define  a  decision  point  for  further  action. 

73.  If  the  Philadelphia  District  considers  restorii^  the  damaged  seaward 
ends  of  the  jetties,  streamlined  jetty  tips  to  reduce  flow  convergence  should  be 
considered. 

74.  Continued  monitoring  should  consist  of  semiannual  hydrogr^hic 
surveys  of  the  bed  and  of  the  jetty  toes  and  bridge  pier  protection  for  2  years, 
then  an  annual  survey  thereafter.  Survey  data  shoidd  be  used  to  update  the 
types  of  data  presentations  used  in  this  report.  Biennial  surveys  of  the  flood 
shoal  will  be  helpful. 

75.  Further  studies  to  prepare  for  contingencies  should  include 
consideration  of  the  following: 

a.  Borings  along  the  inlet  centerline  and  across  Range  7+00  E  to 
define  sediment  strata  down  to  about  100  ft.  Existing  geophysical 
records  should  be  reviewed  for  adequacy  and  supplementing,  if 
needed. 

b.  Erosion  testing  of  the  bed  clays  at  several  dq>ths  to  determine  the 
critical  shear  stress  for  erosion  and  the  erosion  rate  constant  under 
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flow  shear  and  flow  plus  sand  abrasion.  Sample  collection  and 
care  must  carefully  conform  to  laboratory  analysis  specifications. 


c.  Measurement  of  three-dimensional  currents  with  a  boat-mounted 
Acoustic  Doppler  Current  Proflling  (ADCP)  m^er  throughout  the 
inlet  to  siq^ly  data  for  further  analysis  and  to  verify  a  model  if 
needed. 

d.  A  high-resolution,  nonhydrostatic,  three-dimensional  numerical 
model  to  map  bed  stresses,  evaluate  maximum  scour  depths,  and 
evaluate  reni^ial  measures.'  A  physical  model  could  also  be 
used  to  define  the  stresses,  but  c^d  not  adequately  rq)roduce 
scour  and  bed  change  effects  on  hydrodynamics. 

e.  Careful  evaluation  of  any  proposed  remedial  measures  for  other 
impacts,  such  as  those  on  shore  erosion,  water  quality,  or  move¬ 
ment  of  scour  holes  to  new  locations. 


Questions  Posed  by  the  Philadelphia  District 

76.  The  questions  posed  by  the  Philadelphia  District  in  paragraph  3  have 
been  iiiq)licitly  addres^  in  the  preceding  sections.  Here  tte  Committee’s 
views  are  explicitly  listed. 

a.  Wtat  factors  caused  the  accelerated  scour  to  begin  during  the  mid- 

1970s? 

Committee  Comment:  It  is  difficult  to  apportion  the  responsibility  among  the 
several  possible  contributing  causes;  however,  it  is  believed  that  the  erosion 
was  inevitable  once  jetties  were  constructed.  The  accelerated  scour  was  most 
probably  caused  by  exposure  of  the  clay  layers  to  flow  convergence  stress 
concentrations,  but  contributing  factors  could  include  removal  of  the  old 
bridge  and  flood  shoal  mining. 

(1)  Does  sand  bypassing  effect  the  scour? 

Committee  Comment:  The  sand  bypassing  is  not  significantly  affecting  the 
present  scour  rate  since  the  transport  capacity  has  already  stripped  the  sand 
bed.  At  a  significantly  larger  inlet  size,  the  transport  capacity  may  fall  to  a 
point  that  the  100,000  cu  yd  per  year  bypassed  would  otherwise  suffice  to 
maintain  a  stable  bed. 

(2)  Does  the  tidal  prism-inlet  area  theory  apply? 


'  An  ai^rendy  appropriate  model,  MAC-3D,  has  recently  been  developed  by  Dr.  Bob 
Bernard  of  die  Hydraulics  Laboratory.  (Personal  communication,  23  May  1993, 
Dr.  Bernard). 
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Committee  Comment:  The  equilibrium  tidal  prism  curve  of  Jarrett  and 
O’Brien  cannot  be  used  to  directly  infer  equilibrium  inlet  size.  In  combination 
with  the  Dean-O’Brien  stability  curve,  it  can  be  used  to  qualitatively  evaluate 
inlet  stability  under  constant  littoral  siq^ly. 

b.  What  course  enaction  is  recommended? 

Committee  Comment:  Regular  monitoring  and  preparation  of  a  contingency 
plan  by  the  path  described  in  paragraphs  68-72. 

(1)  Is  monitoring  alone  sufficient  at  present  or  should  planning- 
design-budgeting  begin? 

Committee  Comment:  Planning  studies  providing  for  the  analyses  listed  in 
paragnq>h  71  and  72  should  begin. 

(2)  What  is  the  threshold  for  beginning  stabilization  ? 

Committee  Comment:  Plans  for  stabilization  should  proceed  if  the  scour 
expands  its  lateral  extent  toward  the  bridge  or  jetty  toes,  or  if  the  studies 
described  above  indicate  that  scour  will  continue. 

(3)  Are  the  architect/engineer  stabilization  study  findings 
appropriate? 

Committee  Comment:  The  Committee  does  not  consider  the  architect/engineer 
st2d)ilization  schemes  as  listed  to  offer  a  reliable  remedy.  In  particular,  single 
weirs  will  not  be  effective.  An  array  of  weirs  at  several  of  the  locations  listed 
might  stabilize  the  inlet,  but  careful  model  evaluation  will  be  required. 

(4)  Is  local  cross-section  control  sifficient? 

Committee  Comment:  No. 

(5)  Should  present  flows  be  accepted  and  the  bed  armored? 

Committee  Comment:  Reasonable  cost  stabilization  works  are  unlikely  to  alter 
the  flows  enough  to  stabilize  the  inlet.  If  stabilization  works  are  needed,  some 
sort  of  armoring  will  be  the  best  ^proach. 

c.  What  modeling  tools  and/or  prototype  data  are  needed  to  design  a 
stabilization  scheme? 

Committee  Comment:  The  approach  listed  earlier  in  paragraphs  71  and  72  is 
recommended. 

(1)  Should  more  hydrographic  surveys  be  scheduled? 
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Coounittee  Comment;  Yes.  Survey  at  least  semiannually  for  2  years,  then 
annually  thereafter. 

(2)  What  models  are  needed  to  evaluate  the  problem  and  design 
stabilization  features  ? 

Committee  Comment:  See  paragraph  72(d). 
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Table  A1 

Imban  Rivor  Inlot  Hydrography 


-  Range  10+00  East 


Daptha.  ft  Baiow  NOVO,  for  h 


May  88 

Dec  78 

Vr74 

Mar  66 

Nov  60 

Mar  60 

Yr  41 

350.0 

37.2 

36.5 

33.9 

25.9 

10.9 

15.5 

17.5 

10.9 

400.0 

39.1 

29.0 

26.5 

23.7 

10.9 

13.3 

19.0 

10.9 

450.0 

14.1 

11.5 

9.9 

9.8 

10.9 

9.4 

17.5 

13.2 

500.0 


Avg  Depth 


Area,  aq  ft 


16,922  16,227 


7.8 

13.4 

11.6 

9.5 

3,913 

6,700 

6,777 

4,745 

Table  A2 
Indian  River 


Diatanca 
Acroaa 
Range  ft 


Inlet  Hydrography  -  Range  5  +  00  E 


Depths,  ft  Below  NQVD,  for  IntBcated  Data 


Nov  92  Mar  88  Dec  78 


Appendix  A  Hydrographic  Data 


A3 


Tabto  A3 

Indian  Rivar  Inlet  Hydrography  -  Range  0+00 


D^ttha.  ft  B«iow  NQVD.  for  IndteatMl  I 


Dm:  78 

Yr74 

Mar  66 

Nov  60 

Mar  60 

Yr41 

50.0 

24.5 

20.6 

25.5 

19.3 

100.0 

18.5 

41.0 

30.6 

30.2 

150.0 

27.9 

46.7 

41.8 

35.9 

200.0 

54.6 

68.4 

46.5 

34.0 

250.0 

65.7 

70.4 

47.5 

39.7 

300.0 

52.1 

65.8 

44.8 

39.9 

350.0 

28.2 

52.0 

35.8 

34.3 

400.0 

26.9 

38.0 

32.2 

33.1 

450.0 

23.8 

14.0 

12.9 

14.5 

500.0 

-6.0 

-6.0 

-6.0 

-6.0 

Avg  Depth 

28.2 

36.8 

27.8 

24.4 

Area,  sq  ft 

14,100 

18,404 

13,890 

12,222 

TabtoA4 

IncNan  River  Inlet  Hydrography  •  Range  5-t-OO  West 


Dfatanca 

Oaptha,  ft  Baiow  NGVD,  tar  bidtaatad  Data 

Acroaa 

Rangaft 

Nov  92 

May  88 

0ac78 

Yr74 

Mar  66 

Nov  60 

Mar  50 

Yr  41 

o 

b 

-6.0 

-6.0 

-6.0 

-6.0 

-6.0 

-6.0 

-6.0 

-6.0 

50.0 

12.2 

8.0 

8.4 

10.4 

9.9 

9.7 

10.2 

3.2 

100.0 

20.8 

14.1 

14.1 

12.1 

13.0 

9.7 

10.3 

6.4 

150.0 

29.7 

17.9 

27.6 

12.4 

11.9 

9.4 

13.0 

6.9 

200.0 

56.6 

51.6 

34.1 

22.0 

11.3 

7.6 

14.4 

6.6 

250.0 

75.3 

70.9 

42.1 

25.0 

12.8 

8.4 

16.9 

9.3 

300.0 

81.8 

78.5 

41.2 

27.4 

18.4 

8.5 

19.2 

11.2 
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TaMaAS 
IncHan  Rivar  Inlat 


DtetwiM 

AcsroM 

RanSdft 

Nov  92 

o 

6 

-6.0 

50.0 

17.6 

100.0 

27.0 

150.0 

26.5 

200.0 

30.4 

250.0 

33.6 

300.0 

32.0 

350.0 

28.6 

400.0 

27.5 

450.0 

31.6 

500.0 

33.9 

550.0 

33.6 

ISI9II 

33.5 

650.0 

34.6 

700.0 

36.4 

750.0 

40.4 

800.0 

20.7 

850.0 

-6.0 

Avg  Depth 

28.4 

Area,  aq  ft 

24,129 

Hydrography  -  Ranga  10+00  West 


OapIlM,  ft  Bdow  NOVO,  for  hMlcaMd 


Dee  78 

Yr74 

Mares 

Nov  60 

Mar  SO 

Yr  41 

25,255  21,058 


17,021  15,172 


14,375  10,763  5,843 
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